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Abstract 

The g-state Potts model has stood at the frontier of research in statistical mechanics for many 
years. In the absence of a closed-form solution, much of the past efforts have focused on locating 
its critical manifold, trajectory in the parameter {q, e J } space where J is the reduced interaction, 
along which the free energy is singular. However, except in isolated cases, antiferromagnetic (AF) 
models with J < have been largely neglected. In this paper we consider the Potts model with AF 
interactions focusing on obtaining its critical manifold in exact and/or closed- form expressions. 

We first re-examine the known critical frontiers in light of AF interactions. For the square lattice 
we confirm the Potts self-dual point to be the sole critical frontier for J > 0. We also locate its 
critical frontier for J < and find it to coincide with a solvability condition observed by Baxter 
in 1982. For the honeycomb lattice we show that the known critical frontier holds for all J, and 
determine its critical q c = ^(3 + \/5) = 2.61803 beyond which there is no transition. For the 
triangular lattice we confirm the known critical frontier to hold only for J > 0. 

More generally we consider the centered-triangle (CT) and Union-Jack (UJ) lattices consisting 
of mixed J and K interactions, and deduce critical manifolds under homogeneity hypotheses. For 
K = the CT lattice is the diced lattice, and we determine its critical manifold for all J and find 
q c = 3.32472. For K = the UJ lattice is the square lattice and from this we deduce both the 
J > and J < critical manifolds and q c = 3. Our theoretical predictions are compared with 
recent numerical results. 
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The Potts model, proposed by Potts exactly sixty years ago []| , has been at the forefront 
of interest throughout the years. The Potts model extends the 2-state Ising model to q > 2 
states (for reviews on relevances of the Potts model see 2, sj). Despite the intense interest, 
however, the general g-state Potts model has remained unsolved. An equivalent and more 
revealing formulation of the Potts model is the random cluster model advanced by Kasteleyn 
and Fortuin 0, Q . In this formulation the partition function is written as a sum of randomly 
connected cluster of lattice sites with edge weights v = e J — 1, where J is the Potts interaction 
with J > (resp. < 0) indicating ferromagnetic (resp. antiferromagnetic (AF)) interactions. 
The number of spin states q then emerges as a cluster weighting factor, thus providing a 
powerful means of continuing the Potts model to non-integral values of q. In the limit 
of q — > 1, for example, the formulation leads to a bond percolation with edge occupation 
probability p — 1 — e~ J . 

Most studies of the Potts model have focused on ferromagnetic interactions, and the 
AF models have been largely neglected. (See however where the AF Potts model is 

discussed in the context of field theory.) For J < 0, the cluster weights which are products 
of edge weights v, can be either positive or negative resulting in a delicate cancelation. 
Since most analyses are based on a positivity assumption of Boltzmann weights, they are 
not applicable to AF models. In the percolation limit, for example, the edge occupying 
probability becomes negative losing its physical meaning. 

In the absence of closed-form solutions, efforts have focused on determining its critical 
manifold, or critical frontier, a trajectory in the {q, e J } parameter space along which the 
free energy becomes singular. But known critical manifold of the Potts model are very few. 
They are limited to the square, honeycomb and triangular lattices given by 



v 
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q (square) (1) 



v 3 — 3qv — q 2 = (honeycomb) (2) 

v 3 + 3v 2 — q = 0. (triangular) (3) 

Here, we examine closely the applicability of ([1]) - ([3]) to AF interactions. 

The critical point (CD) for the square lattice was obtained by Potts [l| using the duality 



relation [2] 

(e J -l)(e J * -l) = q (4) 
relating an interaction J to a J* on the dual lattice. Since the square lattice is self-dual 



and assuming a unique transition, the duality (jlj) determines the critical point at J = J* or 
v 2 = q. This works only for J > 0. For J < and q > 1, @ relates J to a dual J* in the 
complex-temperature plane. The self-dual argument fails and (DQ) does not hold. 

In 1982, Baxter |8|] observed that the AF Potts model on the square lattice is exactly 
solvable at 

(e J + l) 2 = 4-g, (g<3) (5) 

and argued that the condition ([5]) should coincide with the critical point as in the case of 
the ferromagnetic model. Indeed, as we shall see below, the condition fl5]) emerges as the 
AF branch of the Potts critical manifold when the square lattice is considered in the larger 
setting of a Union- Jack lattice. 

One important feature of the J < transition is the existence of a critical q c beyond 
which there is no transition. The disappearance of a transition for large q in any AF model 
is expected on general ground since for q sufficiently large the ground state entropy will 
destroy any transition. But the actual determination of q c is lattice-dependent. For the 
square lattice, for example, ()5]) gives q c = 3. 

The critical points ([2]) and (|3]) for the honeycomb and triangular lattices can be derived 
as follows: Consider a triangular Potts lattice whose Hamiltonian is separable into sum over 
individual up-pointing triangles so that the partition function can be written as 

Z = Y,X\W{s l ,s 2 ,s z ), (6) 

Si=l A 

where W(s\, s 2 , s%) is the Boltzmann weight of an individual triangle with 3 corner spins 
sij s 2j S3; the product is taken over all up-pointing triangles. An example is the honeycomb 
lattice divided into triangles shown in Figfjja). Clearly, properties of this Potts model are 
completely determined by the weights W(s -t , s 2 , 



Baxter et al 



9] (see also Wu and Lin lOf) have obtained a self-dual trajectory for a 



triangular Potts model with 2- and 3-site interactions. In terms of the ^-weights the self- 
dual trajectory can be recast as 

W(l, 1,1) = 3 W(l, 1, 2) + (q - 2) W{1, 2, 3). (7) 

Wu and Zia [ill have further established that, for positive W(si, s 2 , s 3 ) and in the regime 
W(l, 1, 1) dominates, or 

W(l, 1, 1) > {W(l, 1, 2), W(l, 2, 3)} > 0, (8) 



there is a transition into a ferromagnetic ordered state on the underlying triangular lattice 
at the self-dual point (J7J). Indeed, for W > the self-dual trajectory ([7]) lies entirely in 
the regime (|S]). The transition is expected to be continuous for q < 4 and of first order 
for q > 4 from universality considerations [2|. In the regime where W(l,2,3) dominates 
and for q = 3, there exists a transition into an antiferromagnetic ordered state which has a 
6-fold degenerate ground state. But the critical manifold for this transition in terms of the 
W- weights is not known. 

We now specialize ([7j) to specific lattices. For the honeycomb Potts lattice we take as 
the individual up-pointing triangles the setup in Fig. Ufa). This gives after tracing over the 
center spin sq, 

W(l, 1, 1) = e 3J + q - 1, W{1, 1, 2) = e 2J + e J + q - 2, W{1, 2, 3) = 3e J + q - 3. (9) 

Substituting Q into ([7j) we obtain the critical point (JSJ). Since > for all (positive and 
negative) J, the critical manifold (j2]) holds for all J and gives the critical q c as the root of 

q 2 — 3q + 1 = 0, (honeycomb) (10) 

or q c = |(3 + >/5) = 2.61803.... 

For latter use we recast the critical manifold (|7j) in another language. The star-triangle 
transformation of Fig. [1] transforms the honeycomb lattice to a triangular lattice with 2-site 
interactions K' and 3-site interactions M in every up-pointing triangle. A little algebra 
shows that the critical manifold (JJJ) can be written as 

e ^'+w_ 3e ^_( g _2) = 0, (11) 

with 

e 3^'+M = e 3J + q-l K , = e 2J + e J + q - 2 

3e J + g-3' 3e J + g - 3 1 ' 

Now turning to the triangular Potts lattice for which each individual triangle consists of 

3 edges of interactions J, we have 

W(lA,l) = e 3J , W(l,l,2) = e J , W(l,2,3) = l. (13) 

For J > the weights (fT3|) are in the regime ([8]) and the substitution of (|T3|) into ((?]) leads 
to the critical manifold ([3]). For J < and g = 3, the weights f lT3|) are in the AF regime. 
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The location of its critical point has been the subject matter of several studies. From low 



temperature series of the partition function Enting and Wu 
transition at the critical point e Jc = 0.204. Adler et al 
from Monte Carlo simulations, and Chang et al [ijj found e J ' 



121 1 determined a first-order 
13| found e Jc = 0.20309 ± 0.00003 
0.203073(20) from transfer- 



matrix computations. Feldmann et al 15| showed this transition to be consistent with a 

n 

complex singularity in the dual honeycomb lattice |16| . 





3 K' 



(a) (b) 



FIG. 1. A star-triangle transformation. The 'star' in (a) is mapped into a triangle in (b) consisting 
of 2-spin interactions K' and 3-spin interaction M . 

We consider next the critical point of the diced lattice, a subject matter of recent studies 



171 . 1 181 ] . The diced lattice is the lattice in Fig. |2]^a) with the heavy edges K removed. 
In 1979, one of us [l^ (see also j^oj]) introduced a homogeneity hypothesis and derived 
a conjectured critical manifold for the diced lattice. The derivation of the homogeneity 
hypothesis is as follows: 

Apply the star-triangle transformation of Fig. [1] to all 3-coordinated sites of the diced 
lattice. This results in a triangular lattice with 2-site interactions 2K' and 3-site interactions 
M in every face with e K ' and e M given by f fl2|) . If the interaction M exists in only the 
up-pointing triangular faces, the exact critical frontier is ( fill with K' — > 2K' . The idea of 
homogeneity hypothesis, which we shall later extend to the Union- Jack lattice, is to postulate 
that when 3-site interactions M and M' are present in alternating faces, the critical frontier 
is (jTTj) modified by replacing M — > M + M' in the exponential. In the case of diced lattice 
where M' = M, this gives the critical manifold CP with K' 2K', M 2M. This is the 
homogeneity approximation (HA). 

The resulting critical manifold can be written in a generic form 

W 2 (l,l,l) = 3H/ 2 (l,l,2) + (g-2) W 2 (l,2,3) (HA) (14) 
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with W- weights given in fl9]). Explicitly, ( IT4|) reads 



v e + 6t> 5 + 12t> 4 + 2qv 3 - 9qv 2 - 6q 2 v - q 3 = 0. (HA - diced lattice) (15) 

The critical manifold for the kagome lattice is the dual of (!T5|) . or 

v 6 + 6v 5 + 9v 4 - 2qv 3 - Ylqv 2 - 6q 2 v - q 3 = 0. (HA - kagome lattice) (16) 

This is the conjecture advanced in In the case of M ^ M' and the up- and down- 

pointing triangular faces have different weights and W v , the generalization of (TT4"|) is to 
replace W 2 by Wa ■ W v . 

The critical manifolds f lT5|) and f fTB"]) are exact for the Ising model q = 2. At g = 2 f lT5|) 
gives the solution cosh J = (a/3 + l)/2 and ([TBI) has only one solution at e J = a/3 + 2\/3 = 
2.5424598. both in agreement with Syozi [21 1. 



At q — 3, ( TT5l) has six roots, (e J 



) diced 



u+1 = 2.598917, 0.121599, -2.255914, -0.462026, 



and -0.00128784 ± 1.742380 i. The first 2 roots give respective J > and J < critical 
points. The other 4 solutions are outside the regime (jSj) and are not singular points, so 
they are discarded. The J > transition point e J = 2.598917 is within 0.002% of the 



high-precision value e J = 2.598755(6) from a preliminary finite-size analysis 22j. But the 
J < transition point 0.121599 differs by about 8% from the value 0.1393650(4) from the 
same finite-size analysis. Earlier numerical estimates include 0.1394 and 0.1380 of Kotecky 



et al 



3 



and Chen et al 



and the value 0.1393(8) of Feldmann et al 



from a complex-temperature singularity of the dual kagome lattice 16] . All these values are 



231 ] computed 



consistent with the precise determination 0.139365 of 
of the HA prediction indicates its inadequacy for AF interactions 
At v = — 1 ( fT5i) determines q c for the diced lattice from 



221 ] . The relatively sizable discrepancy 



6g 2 + llg-7 = 0, 



(diced) 



(17) 



or q c = 3.32472. This is compared to the bound q c = 3.117689 argued by Kotecky et al [17] 
and the more precise value q c = 3.4380(2) from finite-size analysis 22]. 

The critical points of the kagome lattice are the dual of those of the diced lattice. For 
q = 3 the singularity at 2.598917 of the diced gives a corresponding kagome critical point 
2.876269 which has been shown 24| to be of the same high accuracy as the diced solution. 



The dual of the AF singularity is complex having no physical meaning. So there is only 



a single (ferromagnetic) transition in the kagome lattice. Note that a direct evaluation of 
(fTEj) yields a spurious AF solution e Jc = 0.078600 which is the dual of the discarded diced 
solution —2.255914 not representing a singularity. 




(a) 




FIG. 2. (a) The centered-triangle lattice constructed by inserting sites (open circles) to faces of 
the triangular lattice, (b) The Union-Jack lattice. 



We consider next the centered-triangular (CT) and Union- Jack (UJ) lattices shown in 
Fig. |2j The CT lattice, also known as the asanoha or hemp- leaf lattice [2l|], is a triangular 
lattice with a spin inserted in each triangular face. Similarly the UJ lattice can be regarded 



as a centered-square lattice. These lattices have received considerab 
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b 



e recent attention as 



25). 



instances of finite-temperature transitions in AF Potts models 

For the CT lattice we trace over all 3-coordinated sites. This leads to a triangular lattice 
with 2-site interactions K + 2K' and 3-site interactions M in every face. Then, as in the case 
of the diced lattice, the homogeneity CT critical manifold is ( ITT]) with K' — > K + 2K', M — > 
2M. This leads to the critical manifold (fT4l) with 



W(1,1A) =e 3K/2 (e 3J + q-l), 
W{1,1,2) = e K/2 (e 2J + e J + q- 2), 
W(l,2,3) =3e J + g-3. 



(18) 



(Note the split of the interaction K into two parallel K/2.) Explicitly, the critical manifold 

reads 



e iK (e- iJ +q-iy = 3e K (e 2J +e J +q-2y + (q-2)(3e J +q-3Y 1 (HA - CT lattice) (19) 

which is to hold in regime (IH|) with W- weights (1181) . 

The critical manifold ( JT9l) is exact for the Ising model since the star-triangle transforma- 
tion is exact for q = 2. This gives 2 cosh J = 1 + \/3 e~ K . In the general case we have 
considered the solution of f )T9|) for \K\ = \ J\. We find q = 3 solutions e J = e K = 1.596063 



and e° 



-K 



0.602542. These values are accurate to within 5 decimal places when 



22|. There is no K < 0, \J\ = \K\ 



compared with results of a preliminary finite-size analysis 
solution. 

The Union- Jack is the lattice shown in Fig. [2](b) . In order to extend the homogeneity 
consideration, we consider first the anisotropic UJ lattice with a unit cell shown in Fig. [3]^a). 
Setting J 4 = as in Fig. Mjo), the UJ lattice becomes the triangular lattice in Fig. H]with 
Ji, J2, J3 interactions in up-pointing triangular faces. The critical manifold of this triangular 
lattice is the anisotropic version of (J7J) [11!], namely. 



W(l, 1, 1) = W(l, 2, 2) + W(2, 1,2) + W{2, 2, 1) + (q - 2) W(l, 2, 3). 



(20) 



Here (120 p is to hold in the regime W > and W(l, 1, 1) dominates, and the transition is to 
a state of ferromagnetic ordering. 

The W- weights in f )20|) can be read off from Fig. [3](b), 



W(l, 1, 1) = e Kl+K2 (e Jl+j2+J:i + q-l), 



W(l, 2, 2) = e K2 (e Jl + e Wz + q - 2), 
W(2, 1, 2) = e Ja + e Jl+Ja + q - 2, 
W (2,2,1) = e Kl (e J:i + e Jl+Ja +q-2), 
W(l, 2, 3) = e Jl + e Ja + e Ja + q - 3, 



(21) 



and the critical manifold (1201) with fl2Tl) is exact. 





FIG. 3. (a) Union-Jack lattice with anisotropic interactions, (b) The J4 = Union-Jack lattice. 

Now consider the Union- Jack lattice with J 4 7^ 0. As in the case of the CT lattice, 
we introduce the homogeneity hypothesis by inserting J 4 into (12 Op such that the resulting 
expression retains appropriate symmetries of the 4 interactions J±, J2, J3, J4, and reduces to 
ff20|) when J 4 is set to zero. The simplest way to do this is to insert J 4 into appropriate 
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FIG. 4. A triangular lattice with J interactions in every up-pointing triangular face. 



exponents in W{1, 1, 1), W{1, 2, 2), W(2, 1, 2), W(2, 2, 1) and insert e h - 1 into W(l,2,3). 
This gives rise to the following homogeneity approximation for the anisotropic UJ lattice, 

e K 1 +K 2 ( e .h+j 2 +j 3 +j 4 + q _ -g = e ^2^Ji+j 4 + e J 2 +j 3 + 5 _ 2) 

+ e Kl (e Js+j4 + e Jl+j2 + q - 2) + e j2+j4 + e Jl+Ja + g - 2 
+ (g - 2)(e Jl + e j2 + e Ja + e j4 + q - 4). (22) 

When Ki = K? = 0, (1221) reduces to the critical manifold of the checkerboard lattice 



conjectured in 19J. However, when K\ ^ K 2 the expression (1221) lacks the intrinsic symmetry 
of J\ -H- J3 or J 2 -H- J 4 . For the isotropic lattice Ji = J2 = J3 = J4 = J we are interested in, 
we obtain 

e ^ 1+ ^( e 4J + g _ 1) = (e*i + e ^)(2e 2J + g - 2) + 2 e 2J + 4(g - 2)e J + (g - 2)(g - 3). 

(HA - Union Jack lattice) (23) 

This is the homogeneity approximation for the Union- Jack lattice. 

The manifolds (1221) and f|2"3"|) are exact at q = 2 as they reproduce the known Ising critical 
point obtained by Wu and Lin [26]. When J = 0, the Union- Jack lattice reduces to a square 
lattice with anisotropic interactions K\ and K2, and (123]) gives the known critical point 
(e Kl — l)(e K2 — 1) = q. When Ki = K 2 = 0, the Union- Jack lattice reduces to a square 
lattice of interactions J and (j2"3"|) factorizes into 

[(e J - l) 2 - q] [e 2J + 2e J + q - 3] = 0. (24) 

Setting the first factor equal to zero gives the critical frontier ([T]). Setting the second factor 
equal to zero gives the Baxter solvability condition (jSJ) as the J < critical frontier as 
alluded to earlier. 

In the general case we have considered the solution of (|23|) for = \J\. For q = 3 we 
found solutions e J = e K = 1.80168 and e J = e~ K = 0.49256. In addition, there is an AF 



solution in 2 < q < 3 at, for example, e K = e J = 0.081504 at q = 2.5. This solution is an 
extension of the J < transition ([5]). For K\ = K 2 = K and setting e K = e J = in (T231) . 
we obtain q c = 3. Indeed, numerically we found no K = J < solution of (1231 in g > 3. 

have recently studied the g = 4, A" = J < 0, Union- Jack 



We remark that Chen et al 



model (see also [25j]) using a tensor-based numerical method, and found an 'entropy-driven' 
transition at e K = e J = 0.0523. Since we found q c = 3 for transitions to ferromagnetic 
ordering on the underlying square lattice, the transition at q = 4 is likely of a different 
nature as described and argued by Deng et al [251 ]. 

Finally we remark that since the duals of the CT and UJ lattices are, respectively, the 3-12 
and 4-8, or (3, 12 2 ) and (4, 8 2 ), lattices, the duals of (IT91) and (1231) give the critical manifolds 
for the 3-12 and 4-8 lattices. In the case of uniform interactions K = K\ = K 2 = J, we 
obtain 

v 7 ( v + 3) 2 = q[3v 7 + (q + 32)v 6 + 15(2g + 5)t» 5 + q(12q + lll)v 4 + 2q 2 (q + Al)v 3 

+ 36gV + 9q 4 v + q% (3 - 12 lattice) (25) 

V Q + 4 V 5 = q ( w 4 + IQ V 3 + I5qv 2 + 6q 2 v + q 3 ), (4-8 lattice) (26) 

where v = e J — 1. Both (1251) and (12B1) are exact at q = 2 reproducing the known Ising critical 



points 2jJ, and have only ferromagnetic solutions. At q — 1 ( 1251) and (|26l) give the bond 



percolation thresholds {2c) 3 _ 12 = 0.740423 and (p c ) 4 _ 8 = 0.676835 in agreement with prior 



numerical estimations 



29] 



In summary, we have considered a lingering unsettled question on critical manifolds of 
Potts models with AF interactions. We have re-examined the known critical frontiers of 
the square, honeycomb and the triangular lattices in light of AF interactions. We have 
also examined the critical manifolds of the diced and kagome lattices under a homogeneity 
assumption, and extended the homogeneity consideration to centered-triangular and Union- 
Jack lattices to deduce critical manifolds. Our theoretical predictions are compared with 
results of several recent numerical studies. 



After the completion of this work, we received a preprint [27J on a new graphical analysis 
of the Potts critical manifold which identifies the kagome conjecture f[T6"|) as a first-order 
approximant, thus resolving the long-standing question on the exactness of the conjecture. 
In 27J the AF critical line (J5J is also deduced from the checkerboard manifold. We thank 



J. L. Jacobsen for sending a copy of [27| prior to publication. The identification of the 
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homogeneity approximation as a first-order approximant in the graphical analysis can also 
be extended to fl23|) for K\ = K2 28j. We would like to thank R. Shrock for comments and 



suggestions and calling attention to 14j-|l6|. and thank J. Salas and R. Ziff for comments. 



We thank Z. Fu for help in the preparation of the manuscript. The work of WG is supported 
by the National Science Foundation of China (NSFC) under Grant No. 11175018, and by 
the program for New Century Excellent Talents in University (NCET-08-0053). 
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